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PREFACE 

This  report  was  written  to  supplement  material  by 
S.  M.  Pollock  in  Selected  Methods  and  Models  in  Military 
Operations  Research  as  well  as  material  in  OEG  Report  No.  56  by 
B.  0.  Koopman.   Much  of  the  report  is  based  on  a  course  given 
by  S .  M.  Pollock  at  the  Naval  Postgraduate  School  in  1969. 
The  report's  relation  to  OEG  Report  56  will  be  evident  on 
reference  to  that  work. 

Some  other  sources  of  material  to  which  this  report 
relates  are  listed  in  the  Bibliography. 
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I.   Detection  Theory  and  Detection  Models 

In  signal  detection  theory,  the  decision  making  portion 
of  a  detection  system  is  called  the  fmczlvzK   and  a  detection 
experiment  is  the  observation  by  a  receiver  of  its  input  during  a 
time  interval.   The  input  which  is  related  to  a  target  is  called 
signal,  and  the  input  which  is  not  related  to  the  target  is  called 
noti><L.       In  general,  the  observation  is  assumed  to  be  of  a  known 
region  which  in  some  cases  is  called  a  /ie-6 o lati.0 n   cull. 

Signal  detection  theory  is  a  basis  for  detection  modeling. 
The  detection  models  which  are  discussed  here  rely  heavily  upon 
it.   When  a  detection  experiment  has  been  performed  either  the 
event  H,  =  {At  least  one  target  was  present  in  the  region  which 
was  observed  during  the  time  of  the  observation.}   or  its  comple- 
ment  H^   will  have  occurred.   In  general,  the  models  will  specify 
that  either  the  event   D,  =  {The  receiver  decides  at  least  one 
target  was  present  in  the  region  which  was  observed  during  the 
time  of  the  observation.}   or  its  complement   Dn   will  have  occurred. 
In  terms  of  signal  and  noise,  the  events   D„   and   D,   can  be 
expressed  as  follows:   D~  =  {The  receiver  decides  its  input  during 
the  time  of  the  observation  was  noise.}   and   D,  =  {The  receiver 
decides  its  input  during  the  time  of  the  observation  was  signal 
and  noise.}   The  theory  which  is  the  basis  for  detection  models 
which  have  the  above  properties  is  called  binary  detection  theory. 

Four  events  which  are  important  in  binary  detection  theory 
are  indicated  in  the  Venn  diagram  of  Figure  1. 


Figure  1.   Four  events  important  to  binary  detection  theory. 

The  Venn  diagram  emphasizes  a  decision  problem  which  is 
associated  with  a  receiver.   The  problem  is  this:   Under  what  con- 
ditions should  the  event   D,   occur,  that  is,  under  what  conditions 
should  the  receiver  decide  that  a  target  was  present  during  the 
time  of  an  observation?   Detection  theory  may  be  able  to  provide  a 
solution  to  this  problem.   If  this  is  the  case,  the  noise  will 
generally  be  described  as  a  random  process  and  the  signal  will  be 
described  as  either  a  deterministic  or  a  random  process. 

The  following  notation  and  terminology  will  be  used: 
pf  =  P(D, |Hq)   will  be  called  the  false  alarm  probability, 
p,  =  P(D1|H1)   will  be  called  the  detection  probability.   P  =  P(H,) 

-11  be  called  the  prior  probability.   P   is  the  probability  that 
a  target  will  be  in  the  region  to  be  observed  during  the  time  of 
the  observation. 

The  input  to  a  receiver  is  assumed  to  be  a  quantity  whose 
square  is  proportional  to  power.   The  input  at  some  time   t.   will 
be  symbolized  by  y(t.).   The  noise  at   t.   will  be  symbolized  by 


n(ti)   and  the  signal  by   s(t.).   The  noise  and  signal  processes 

will  be  represented  by   {n(t.),t.et}   and   {s(t.),t.et}   where   t 

represents  the  time  interval  during  which  the  receiver 

observes  the  input.   Often   s(t.)   and   n(t.)   can  be  assumed  to 

be  additive  in  which  case   H0  =  {y (t . )  =  n (t . ) , t . £t }   and 

Hl  =  (y(ti)  =n(t±)  +  s(t±)  /tiet>. 

Because  of  the  finite  quantity  of  information  present  at 

the  input  of  a  receiver,   y   needs  to  be  measured  at  only  a  finite 

number  of  points  in  time  in  order  to  be  adequately  determined  over 

an  observation  interval.   For  this  reason,  the  noise  and  signal 

processes  can  be  represented  by   (n(t.)/   i  =  l,...,m)   and 

{s(t.),   i  =  l,...,m}   where   t, ,  ...,t   are  in  the  observation 
1  1      m 

interval  of  length   t.   The  noise  process  is  then  defined  by  a  set 
of  m   random  variables.   To  specify  the  noise  process,  one  needs 
to  specify  only  the  joint  distribution  of  the   m   random  variables, 
If  the  signal  process  is  deterministic,  it  is  simply  a  finite  set 
of  values  which  is  known  before  the  detection  experiment  is  per- 
formed. 


II.   Decision  Criteria 

To  simplify  the  discussion  of  decision  criteria  and  decision 
rules,  the  stochastic  process  which  represents  a  receiver's  input 
will  be  assumed  initially  to  be  a  single  random  variable   Y.   It 
can  be  called  the  decision  random  variable.   The  input  process  in 
this  case  is  determined  by  the  two  conditional  distribution  func- 
tions  F  (y|HQ)   and   F  (y|H,). 

The  condition  that  a  receiver's  input  is  required  to  satisfy 
in  order  that   D,   will  occur  can  be  specified  in  terms  of  a  deci- 
sion rule.   For  the  assumed  case,  a  decision  rule  is  a  rule  which 
determines  for  every  observable  value  of   Y   the  decision  that  the 
receiver  is  to  make.   The  decision  rule  can  be  considered  to  be  a 
function   <p  (y)   which  relates  the  observable  values   y   to  the 
following  two  statements: 

d~:   Decide  the  input  was  noise, 
d, :   Decide  the  input  was  signal  and  noise. 
Defining  a  decision  rule   ip  (y)   is  equivalent  to  defining  a  set 
tt      such  that   D,  =  {YeQ}. 

The  problem  which  was  considered  above  can  now  be  restated 
in  the  following  way:   What  criterion  should  be  adopted  in  order 
to  determine  a  decision  rule?   A  desirable  characteristic  for  a 
criterion  is  suggested  by  the  following  argument:   Consider  the 
ratio  of  posterior  probabilities 


P(H1|y)/P(H0|y) 


One  might  expect  that  the  values  of   y   in   ft   should  make  this 
ratio  relatively  large.   But  making  this  ratio  large  is  equiva- 
lent to  making  the  likelihood  ratio   L(y)   large.   This  suggests 
that   ft   might  be  defined  as  follows:   ft  =  {y:  L(y)  £  K>   where 
K   is  some  positive  constant  which  has  yet  to  be  specified. 

Four  specific  decision  criteria  are  defined  next  in  terms 
of   ft.   For  each  criterion,  ft   has  the  above  form  with  only  the 
procedure  for  determining  K   being  different.   The  decision  cri- 
teria are: 

1.  The  Neyman-Pearson  Criterion:   Choose   ft   so  that   p,   is 
maximized  subject  to  the  constraint  that   pf  £  a   where   a   is 

a  specified  value.   For  a  continuous  decision  random  variable 
the  constant   K   is  chosen  so  that   p.  =  a. 

2.  The  Bayes  Criterion:   Choose   ft   so  that  the  expected  cost 
of  a  receiver's  decision  is  a  minimum.   For  a  continuous  decision 
random  variable   K  =  (C10-C00/C01-C1;L)  (1-P)/P   if   ci0  >  C00   and 

C-.  >  C, ,   where   C.    is  the  cost  of  the  event   D.  fi H . . 
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3.  The  Ideal  Observer  Criterion:   Choose   ft   so  that  the 
probability  that  the  receiver  makes  an  incorrect  decision  is  a 
minimum.   For  a  continuous  decision  random  variable   K  =  (1-P)/P. 

4.  The  Minimax  Criterion:   Choose   ft   so  that  the  maximum 
expected  cost  of  a  receiver's  decision  is  a  minimum.   For  a  con- 
tinuous decision  random  variable,   K  =  (C,  o~coo/'C01~<~ll^'''~P*^ //p* 
if   C  Q  >  CQ0   and  CQ,  >  C,,.   Here   P*   is  the  value  of   P   which 
would  make  the  expected  cost  of  a  receiver's  decision  a  maximum 

if  a  Bayes  decision  rule   cp    were  used. 


A  more  general  discussion  of  the  above  criteria  would  involve 
the  notion  of  a  randomizing  rule. 

If  a  model  is  adopted  which  specifies  the  conditional  dis- 
tributions  and  a  decision  rule,  then  the  values  of   pf   and   p. 
can  in  principle  be  determined.   The  pair  of  values   (pf/p, )   is 
called  a  receiver  operating  point.   If  the  decision  rule  results 
from  one  of  the  four  likelihood  ratio  criterion  listed  above,  then 
it  will  involve  the  parameter   K   through  the  relation 

ft  =  {y:  L(y)  £K}.   And,  for  a  given  value  of   K,   since   ft   determines 
(PfrPfl),    a  unique  operating  point  result.   By  varying   K,   a  set 
of  operating  points  can  be  generated  which  is  called  an  ROC  curve. 
Different  ROC  curves  can  be  produced  by  changing  one  or  both  of 
the  conditional  distributions.   Note,  a  change  in  a  conditional 
distribution  function  implies  a  change  in  the  signal  or  the  noise. 

Decision  rules  which  result  from  applying  the  four  likelihood 
ratio  criteria  in  a  model  in  which  the  input  process  is  determined 
by  a  set  of  m  random  variables  can  be  expressed  in  terms  of  a  set  ft 
as  follows :   ft  =  { (y_ , . . . ,y  ) :   L(y, , . . . ,y  )  £  K}   where   K   is 
specified  in  the  same  way  as  it  is  in  the  corresponding  case  in 
which  the  input  process  is  determined  by  a  single  random  variable. 


III.   Two  Detection  Models 

A  particular  detection  model  will  now  be  examined.   In 
the  model,  the  receiver's  input  process  is  defined  by: 

H0:   {y (t . )  =  n (t . ) /   i  =  l,...,m}   and 

H^   (y(ti)  =  n(ti)  +  s(t±)  ,   i  =  l,...,m) 

where  the  values  of   y   are  measured  every   At   units  of  time. 

The  noise  process   {n(t.)>   i  =  l,...,m}   is  assumed  to  consist 

of  a  set  of  m   independent  normal  random  variables  each  with 

mean  zero  and  variance   a2.   The  signal  process   {s(t.),  i  =  l,...,m} 

is  assumed  to  be  deterministic.   Thus,  the  input  process  consists 

of  m   independent  normal  random  variables   Y, , . . . ,Y   each  with 

variance   a2   and  each  with  mean  zero  when  the  target  is  not 

present  and  with  mean   s.  =  s(t.)   when  the  target  is  present.   Note, 

such  a  model  might  be  used  to  obtain  an  optimistic  estimate  of  a 

detection  systems  performance,  since  all  the  information  about 

the  signal  is  assumed  to  be  known. 

The  result  of  the  application  of  a  likelihood  ratio  decision 
rule  in  the  model  can  be  expressed  in  terms  of  a  random  variable 
Z.   This  random  variable  is  called  the  cross  correlation  statistic 

and  it  is  defined  by   Z  =  Y  s.Y..   However,  it  is  more  convenient 

u      1  1 

to  express  the  results  in  terms  of  a  random  variable  V  which  is 
defined  by  V  =  Z/a  .   In  terms  of  this  random  variable, 


and 


pf  =  1  -  <D(v*) 


Pd  =  1  -  •Mv*-^) 


where   $   symbolizes  the  standard  normal  distribution  function  and 
v*  =  (l/o)  (a2  lnK+  (1/2)  I   s .  )  ,  ah    =   I    s2/a  . 

The  parameter  d   is  called  the  detection  index  and,  since 
a2  =  I   s2a2,   d  =  I  s2/a2. 

Recall,  the  input  stochastic  process  is  assumed  to  repre- 
sent a  quantity  whose  square  is  proportional  to  the  input  power. 
Therefore,  the  average  power  input  to  the  receiver  over  the  time 
interval  t  during  which  the  receiver  observes  the  input  from  a 
region  can  be  approximated  by  £  y2/m.  The  average  signal  power 
can  be  approximated  by  S  =  £  s2/m  and  the  average  noise  power 
by  I  n2/m.  And  the  expected  average  noise  power  can  be  approxi- 
mated by   N  =  I    a2/m  =  a2 .   in  these  terms,   d  =  m(S/N). 

The  receiver's  bandwidth  will  be  represented  by   BW. 
Assuming  m  corresponds  to  the  value  determined  by  the  sampling 
theorem,  that  is,   m  =  t/At  =  2t(BW) ,   the  detection  index  can  be 
written  as   d  =  2t(BW)(S/N)  or,  defining   NQ  =  N/Bw,   as   d  =  2t(S/NQ) 

A  model  in  which  the  signal  process  is  not  deterministic 
will  be  considered  next.   The  model  is  defined  by 

HQ:   {yi  =  nif   i  =  1, . . . ,m} 

H1:   {yi  =  ni  +  si,   i  =  l,...,m}. 

Here  the  noise  process  is  a  set  of  m   independent  normal  random 

variables  with  mean  zero  and  variance   a2   as  before.   However, 

now  the  signal  process  is  also  a  set  of  m   independent  normal 

random  variables  with  mean  zero  but  with  variance   a2.   The  result 

s 

of  applying  a  likelihood  ratio  decision  rule  in  this  model  can  be 


expressed  in  terms  of  a  random  variable   X  which  is  defined  by 

X  =  £  Y?.   When  a  target  is  not  present,  the  statistic   X/a2   has 

a  chi-square  distribution  with  m  degrees  of  freedom  and  when  a 

target  is  present,  the  statistic   X/(a2+a2)   has  a  chi-square 

distribution  with  m   degrees  of  freedom.   Hence,   pf  =  P(x2^x*/a2) 

and   p,  =  P [x2^x*/ (o2+oz ) }   where   y2   is  a  chi-square  random 
d      m         s  m 

variable  with  m  degrees  of  freedom  and   x*   is  a  number  which 

is  determined  by  the  decision  rule. 

When  the  target  is  not  present,  since   X/a2   has  a  chi-square 

distribution,  the  variance  of   X   is  2moh      and  the  mean  is   ma2. 

When  the  target  is  present,  since   X/(a2+a2)   has  a  chi-square 

distribution,  the  variance  of   X   is   2m(a2+a2)2   and  the  mean  is 

s 

m(a2+a2).   By  the  central  limit  theorem,  as  the  number  of  degrees 
of  freedom  m  of  a  chi-square  random  variable  becomes  large  it 
can  be  approximated  by  a  normal  random  variable  with  the  same  mean 
and  the  same  variance.   Hence,  for  sufficiently  large   m,   pf 
and   p,   can  be  approximated  as  follows: 

Pf  s  1  -  $[(x*-mo2)/(2mo")li] 
and 

Pd  a  1  -  <H  [(x*-ma2)  -ma2]/(2ma,,)^] 

where  (o2+o2)2      has  been  approximated  by   a1*.   This  approximation 
can  also  be  written  as: 

pf  =  1  -  $(v*) 

pd  ~  l  -  $(v*-d^) 
with  v*  =  (x*-ma2)/(2ma'+)^    and   d  =  (ma2)  2/2ma1* . 


Note,  over  the  observation  interval   t,   £  s2 /m   is  approximately 
the  average  signal  power  so   S  =  £  a2/m  =  a2   is  approximately 
the  expected  average  signal  power  and,  as  before,   N  =  a2   is 
approximately  the  expected  average  noise  power.   Hence,  for  this 
model  the  detection  index  can  be  written  as 

d  =  j  (S/N)2  =  t(BW)(S/N)2. 
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IV.   Detection  Model  Applications 

When  a  likelihood  ratio  decision  rule  was  used  with  either 
of  the  two  models  discussed  above,  the  following  result  was  obtained 

pf  =  1-  *(v*) 

Pd  =  1  -  $(v*-d  ) 

In  the  case  of  a  known  signal  (deterministic  process) ,  a  cross 
correlation  receiver  is  required,  and  d  =  2t(BW)(S/N).   This  will 
be  called  Case  I.   In  the  case  of  a  Gaussian  signal,  a  square  law 
receiver  is  required  and,  for   S/N  <<  1   and   (BW)  large, 
d =  t(BW) (S/N) 2 .   This  will  be  called  Case  II.   In  both  Case  I 
and  Case  II,  the  detection  index   d   is  a  function  of   (S/N). 
The  two  models  can  be  related  to  radar  and  sonar  systems,  for 
example,  by  using  the  radar  equation  or  sonar  equation  to  deter- 
mine  (S/N) . 

In  some  radar  and  sonar  detection  models,  it  is  assumed 
that  the  decision  rule  used  results  in  a  required   p_   and  in 
addition  that  a  minimum  value  of   p,   is  specified  in  the  sense 
that  a  target  is  said  to  be  detectable  only  when   p,   is  greater 
than  or  equal  to  the  minimum  value.   This  minimum  value  of   p, 
along  with  the  required   p   define  what  could  be  called  a  mini- 
mum acceptable  signal-to-noise  ratio   (S/N)m-   In  sonar  models,  this 
minimum  acceptable  signal-to-noise  ratio  determines  the  detection 
threshold   DT   by  the  relation   DT  =  10  log (S/N)  .   If  the  mini- 
mum acceptable  value  of   p,   is  .5,  then   DT   is  usually  called 
the  recognition  differential   RD. 
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A  passive  sonar  detection  model  will  now  be  discussed.   It 
an  example  of  a  model  which  uses  the  concept  of  a  minimum  value 
p, .   It  also  illustrates  a  method  of  dealing  with  non-station- 
ry  noise  and  signal  processes.   In  the  model,  the  event  detection 

the  event   {X   £  0}   where   X    represents  signal  excess,  a 
andom  variable  with  expected  value   SE.   The  source  level,  noise 
evel,  directivity  index,  and  the  detection  threshold  or  recogni- 
i-on  differential  are  also  random  variables  with  expected  values 
NL,   DI ,   and   DT   respectively.   The  passive  sonar  equation 
be  written  in  terms  of  these  random  variables  or  in  terms  of 
eir  expected  values.   In  the  latter  terms,   SE  =  FOM  -  TL,   where 

SL  -  (NL-DI)  -  DT   is  the  figure  of  merit.   To  determine  the 
robability  of  a  detection  P  (X„„2>0)  ,      the  distribution  of  the 
adorn  variable   Xor,  must  be  specified. 

This  model  could  be  considered  to  be  equivalent  to  one  in 

ich  the  event   {X   ^0}   is  the  event  that  p,   is  greater  than 

minimum  acceptable  value.   Although  this  definition  may  appear 

Dre  reasonable  than  the  one  which  defines   {Xor,^0}   as  the  event 

Section,  it  would  likely  be  no  more  useful  for  developing  tac- 

s,  designing  exercises,  or  testing  detection  devices. 

The  noise  and  signal  processes  should  be  effectively  sta- 
onary  during  a  "look"  at  a  resolution  cell,  that  is,  over  an 
tegration  period  if  a  Case  I  or  Case  II  model  is  used  to  relate 
p, ,   and  the  signal-to-noise  ratio.   If  either  of  these  models 
used,  the  signal-to-noise  ratio  for  an  integration  period  will 
the  value  of  a  random  variable,  and,  therefore,   p,   will  be 


12 


conditioned  on  the  random  variable  signal-to-noise  ratio.   The 
receiver  is  usually  assumed  to  be  readjusted  so  that   pf   remains 
relatively  constant  from  "look"  to  "look"  at  a  given  resolution 
cell  in  applications  of  the  model. 
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V.   Search  Theory  and  Search  Models 

Search  theory  deals  with  problems  such  as  that  of  determining 
the  optimal  allocation  of  search  effort  in  a  region  or  that  of 
determining  the  probability  of  detecting  a  target  within  a  given 
length  of  time  for  a  particular  search  plan.   In  treating  these 
problems,  false  alarms  are  often  ignored.   In  effect,  the  time 
spent  on  false  alarms  is  assumed  to  be  small  relative  to  the  time 
scale  of  the  problem,  however,  the  cost  associated  with  a  false 
alarm  is  not  assumed  to  be  negligible. 

In  general,  it  is  assumed  here  that  during  a  search  a 
target  will  move  through  resolution  cells  surveyed  by  the  detec- 
tion system  of  a  searcher.   In  the  model  of  such  a  search,  the 
probability  that  the  input  from  a  resolution  cell  containing  the 
target  will  cause  the  receiver  to  declare  the  target  is  present 
is   p , .   In  systems  with  operators,  the  operator  can  delay  the 
decision  regarding  the  presence  or  absence  of  a  target  in  a  reso- 
lution cell  and  can  recall  the  input  of  adjacent  resolution  cells. 
Such  detection  systems  are  difficult  to  model.   To  illustrate 
this,  consider  the  following  hypothetical  system  consisting  of  an 
active  sonar  system  plus  an  operator.   Assume  that  p,   as  a  func- 
tion of  the  signal-to-noise  ratio  has  been  determined  for  the 
operator  in  a  laboratory  experiment  by  forcing  the  operator  to 
respond  after  a  single  "look"  at  a  resolution  cell.   Operationally, 
the  probability  that  the  operator  detects  a  target  will  not  be 
given  by  the  laboratory  determined  value  of   p,   if  the  operator 
delays  the  detection  decision  for  several  "looks."   Various  models 
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have  been  proposed  to  deal  with  this  situation.   In  one,  a  "three- 
out-five"  detection  criterion  has  been  adopted.   The  criterion 
states  that  an  operator  will  declare  a  target  if,  out  of  five 
consecutive  inputs  from  adjacent  resolution  cells,  at  least  three 
are  such  that  they  would  have  caused  the  operator  to  declare  a 
target  present  in  a  single  "look"  experiment.   An  input  which 
satisfies  this  condition  is  referred  to  as  a  success.   With  this 
criterion,  the  probability  that  a  target  will  be  detected  with 
m   consectuive  "looks"  at  the  target  is  found  as  follows:   Deter- 
mine the   2    sequences  of  successes  and  failures  of  length  m 
which  could  result  with  m  consecutive  "looks."   The  probability 
is  equal  to  the  sum  of  the  probilities  of  occurrence  of  those 
sequences  for  which  the  three-out-of-f ive  criterion  is  satisfied. 
The  probability  of  a  particular  sequence  will  depend  on  the  value 
of   p,   and   1  -  p ,   for  each  of  the  m   "looks"  in  the  sequence. 
Unfortunately,  in  most  situations,  the  number  of  sequences  to  be 
considered  is  too  large  to  be  computationally  tractable  even  if 
the  probabilities  in  the  sequence  can  be  determined. 
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VI.   The  Probability  of  Detection  During  a  Search 

The  problem  of  determining  the  probability  that  a  target 

will  be  detected  in  a  search  at  or  before   n   "looks"  have  occurred 

is  generally  basic  to  the  solution  of  search  problems.   Let  N 

represent  the  number  of  the  "look"  at  which  detection  first  occurs, 

then  this  probability  can  be  expressed  as 

n 
P(N^n)  =  I     P(N=i) . 
i=l 

It  can  also  be  expressed  as 

n 

p(N*n)  =1-77  (i-g,) 

i=l 


where  g.  =  P (N=i | N^i-1)  is  the  probability  that  the  event  detec- 
tion  occurs  on  the  i —  look  conditioned  on  the  event  detection 
has  not  occurred  earlier.  The  second  expression  for  P(N^n)  is 
generally  of  greater  interest  than  the  first  one,  since  g.  can 
be  directly  related  to  operational  parameters  such  as  target  range. 
Note,  if  the  resolution  cell  being  examined  on  the  i —  look  is 
empty,  then   g.   will  be  zero. 

A  continuous  analog  to   P(N^n)   can  be  developed  as  follows: 
Consider  a  model  in  which  the  time  for  a  look  is   At,   and  if 
detection  occurs  on  the   i —   "look,"  the  time  of  detection  is 
iAt.   The  probability  of  detection  at  or  before  the   n —  "look"  then 
can  be  expressed  in  terms  of  time  by  noting  that   P(N^n)  =  P(T^nAt) 
where   T   represents  the  time  of  detection.   Now  define  the  detec- 


tion rate  to  be   y (nAt)  =  (1/At) P [T=nAt | T^ (n-1) At]\   The  probability 
P(T^nAt)   can  be  approximated  in  terms  of  a  related  quantity  y(t)   by 
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ft 

P(T^t)  =  l-exp{-     Y(T)dT> 

0 

where   Y(t)   is  related  to  y  (nAt)   by  the  requirement  that 
P(T^t)   equals   P (N*n)   for   t  =  nAt.   In  some  models  where  it 
would  be  useful  to  use   P(T^t),   the  continuous  approximation  to 
P(N£n),  it  might  not  be  reasonable  to  model  the  detection  process 
in  terms  of  single  "look"  probabilities.   If  this  is  the  case, 
Y (t)   must  be  defined  by  some  other  means.   The  visual  detection 
model  developed  in  the  OEG  Report  No.  5  6  is  an  example  of  this. 

Generally,  in  search  operations,  the  target  will  be  moving 
on  a  relative  track   C.   In  some  cases,  it  may  be  useful  to  con- 
sider the  track  to  consist  of  segments   C, ,C2,...,C  .   Let   I*. 
represent  the  event  the  target  is  first  detected  on  track  segment 
C .   The  probability  of  first  detecting  a  target  while  it  moves 
along   C   can  then  be  expressed  as  follows: 

n 
P(T)  =  l-exp{-   I   F(C.)> 

i=l 

where   F(C.)   is  called  the  sighting  potential  on  the  track  seg- 
ment C  .   and  exp{-F(C  )  }  =  P  (F .  |  f .  ,  H  . .  .  fl  f^)  .   It  is  conven- 
ient to  use  the  concept  of  sighting  potential   in  the  analysis  of 
some  search  types,  for  example,  "ladder  searches". 

The  track  of  a  target  which  moves  in  a  plane  can  be  described 
by  equations   x  =  x(t)   and   y  =  y(t).   These  equations  describe 
the  relative  track  of  the  target  with  respect  to  the  searcher 
(often  called  the  observer)  if   x   and  y   refer  to  a  coordinate 
system  in  which  the  searcher  is  located  at  the  origin.   In  this 
case,   r   the  range  of  the  target  is  given  by   r  =  [x2+y2]  2. 
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If  in  a   model  the  detection  capability  of  a  searcher 
against  a  target  is  assumed  to  be  independent  of  the  time  at  which 
the  target  occupies  a  particular  position  on  the  target's  track, 

then  time  can  be  replaced  by  target  distance  along  the  track. 

ft 
This  can  be  done  by  using  the  transformation  equation   s  =    w(x)dT 

J0 
where   w   is  the  target's  speed  relative  to  the  searcher  and  the 

zero  of  time  has  been  chosen  to  coincide  with  the  start  of  the 

search.   If,  in  addition,  the  detection  capability  of  the  searcher 

against  the  target  is  assumed  to  depend  only  on  the  range  of  the 

target,  then  the  sighting  potential  for  a  target's  track  C   for 

a  continuous  model  can  be  expressed  in  terms  of  either   t   or   s 


as  follows:   F(C)  = 
dt/ds  =  l/w(s) . 


ft 

Y[r(T)  ]dT  = 
0 


s 

Y [r (a) ]da/w(a)   since 
0 
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VII.   Experimental  Validation  of  Detection  Models 

Consider  a  detection  system  which  could  be  satisfactorily 
modeled  by  a  continuous  detection  rate  function  that  depends  only 
on   r,   the  target's  range,  if  a  suitable  detection  rate  function 
could  be  found.   Suppose  on  some  basis  that  a  particular  detec- 
tion rate  function  is  proposed?   An  estimate  of  the  suitability 
of  the  proposed  function  could  be  made  using  standard  statistical 
tests  if  appropriate  operational  data  could  be  obtained.   For 
example,  the  number  of  detections  at  various  ranges  for  straight 
line  encounters.   If  such  data  were  available,  its  agreement  with 
the  cumulative  distribution  function  of  the  target's  range  at 
detection  deduced  from  the  proposed  detection  rate  function  could 
be  determined  by  use  of  the  Kolmogorov-Smirnov  test.   This  technique 
could  be  applied  even  though  the  number  of  encounters  in  which  a 
target  is  not  detected  was  unknown.   In  this  case,  the  cumulative 
distribution  function  to  be  compared  is   F  (r|D)   where   D  = 
{detection}   and   R  represents  the  target's  range  at  detection. 
To  derive   F  (r|D)   or  equivalently  the  density  function 
f  (r|D)   given  a  continuous  detection  rate  function,  one  can  pro- 
ceed  as  follows:   First  find  the  joint  density  function  which  de- 
scribes the  distribution  of  the  target's  rectangular  coordinates 
(X,Y)   at  its  detection  point.   To  do  this,  set   f  y(x,y) |D)AxAy  = 
f    (x,y)AxAy/p  where   p  =  P(D).   This  can  be  done,  since 

A,  X 

[  (Y=y)  fid]  =  (Y=y)  .   That  is,   Y   can  be  considered  to  take  on  a 
value  only  when  a  detection  occurs.   Note,   p   is  the  value  of  the 
integral  of   fv   (x,y)   over  all  pairs  of  values   (x,y) .   Now  note  that 

A,  Y 
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fv  „(x,y)    =  fv I v (y I x) fY (x)   and  that  the  conditional  density 

X ,  Y  x  I  X  '       X 

f  I  (y|x)   can  be  written  as 

fYjx(y|x)  =  (l/w)Y[r(y)]exp{-j  y [r (n) ] dn/w} . 

—00 

Given   fv(x)   and   Y(r),   f  (r|D)   can  be  found  by  first  trans- 

X  K 

forming  to  r  and  a  by  using  the  transformation  equations  x  =  r 
cos  a  and  y  =  r  sin  a  and  then  finding  the  marginal  distribution 
of   R   as  follows: 

fR  E  (r,a)da. 


fR(r)  = 


0 
An  example  of  the  above  procedure  will  now  be  given.   Suppose 

l/2a    |x|  £  a  /  kh/r3    |x|  £  a 


=  j 

Ixl  >  a  (    0      |x I  >  a 


fx(x)  =  <  and    y(r) 


The  transformation  to   r   and   a   gives 

r  sin  a 


(l/2awp) (kh/r2)exp{- 

0         |r  cos  a  I  >  a 


kh(r2cos2+n2)  3/2dn/w} 


fR/E(r/a|D)=  }  Ircosal  ^  a 


For  large   a,   this  can  be  simplified  by  letting   a  -*■  °°,   after 

f  a 
noting  that  p  =  l/2a   p(x)dx  =  W/2a   where   W   is  the  sweep 

width.   In  the  limit,  "a  fD  „(r,a|D)  =  (kh/wWr2 ) expl-kh/wr^l-sin  o) } , 

and  after  performing  the  integration  over  a, 

fD(r|D)  =  (W/2r2) (l-$(kh/wr2)^}. 

Note,   fD(r)  =  pf_,(r|D)   which  emphasizes  the  fact  that  in 

general  detection  need  not  take  place  during  a  straight  line  encounter, 
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VIII.   Sweep  Width  Determination  for  a  Random  Track  Angle 

Now  consider  the  problem  of  determining  probability   p  of 
detecting  a  target  during  an  encounter  which  could  be  described 
by  a  straight  line  encounter  model  if  the  angle   cp   between  the 
target's  track  and  searcher's  track  were  known.   Suppose  it  can 
be  assumed  that  lateral  range  and  the  relative  speed 

9      9  *o 

w  =  [u  +v  -2uvcos  cp]  are  values  of  random  variables  whose  joint 
distribution  is  known.  The  probability  p  of  detecting  a  target 
during  a  straight  line  encounter  can  then  be  expressed  as 


.00    .00 


p  = 


p(x,w)fx  ^(x,w)dxdw 


_oo  0 


where   p(x,w)  =  p  (det  |x  =  x  H  ¥  =  w)   and   f   „,(x,w)   is  the  joint 
distribution  of  the  lateral  range   X  and  the  relative  speed   ¥. 

If   X  and   ¥  are  independent  random  variables  and  if 
p(x,w)   as  a  function  of  w   is  adequately  described  by  a  linear 
approximation,  then   p   can  be  approximated  as  follows: 


P  = 


p  (x,w)  fx(x)dx 


where 


w  =  E(iJ0  = 


0 


wfy(w)dw  = 


2tt 


0 


(u2+v2-2uv  cos  (p)  2  f$((p)d(p. 


An  average  sweep  width  can  be  defined  by 

.00 

W  =    p(x,w)dx. 

■ 
—  00 

If   $  is  uniformly  distributed  so   f$(<p)  =  Yt\  '      -  is  the  value 
of  an  elliptic  integral  of  the  second  kind.   An  average  sweep  width 
based  on  this  value  of   w   should  be  used  in  the  Random  Search 
Formula  when  either   u   or  v   are  not  appropriate. 
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IX.   A  Parallel  Sweep  Search  Model 

A  model  of  a  parallel  sweep  or  "ladder"  search  will  be 
considered  next.   In  the  model,  the  target's  velocity  is  zero 
and  its  position  is  uniformly  distributed  in  a  rectangular  region. 
The  searcher  moves  along  a  series  of   n   parallel  tracks  to  that 
n   straight  line  encounters  will  occur  in  a  complete  search  of 
the  region.   The  track  spacing  is   s   and  the  tracks  are  numbered 
from  1  through   n.   The  search  geometry  is  shown  in  Figure  2. 


n 


— r 
ns 


->  z 


Figure  2.   Search  geometry  for  a  parallel  sweep  model. 


If  the  target's  horizontal  coordinate  is  a  random  variable 
Z   in  the  coordinate  system  shown  in  the  figure,  then   P(det|z=z) 


n 


=  l-exp[-      F(x.)]   where   x.  =  z-  (j~h)s  is  the  target' 


j=l 


.th 


lateral  range  from  the   j —   track  and 


/•ns 


P(det)  = 


P(det | Z=z)dz/ns. 
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In  some  cases,  it  is  convenient  to  use  the  following  coordinate 

system:   Let   I   be  the  random  variable  which  is  the  number  of 

the  track  closest  to  the  target.   Relabel  the  tracks  so  that  the 

number  of  the  track  closest  to  the  target  is  0,   then  x.  =  -js  +  x 

where  now   j  =  -K, . . . ,  0, . . . ,M  where   n  =  M  +  K  +  1   and  the  value 

of   K   and  M   are  determined  by  the  value  of   I,   that  is,  by 

the  location  of  the  target.   Now 

m 
P(det|x  =  x  ni=i)  =  l-exp[-   £   F(x.)]  =  p(x,i) 

j=-k    D 


and 


t 
P(det|l  =  i)  = 


s/2 

p(x,i)dx/s  =  p(i) . 

-s/2 


Hence, 


n  n 

P(det)  =  I      p(i)P(I=i)  =  (1/n)  I      p(i) 
i=l  i=l 


or 


n   rs/2 

P(det)  =  (1/ns)  I  p(x,i)dx 

i=l  J  /0 
-s/2 

If   n   is  large  and  the  detection  law  is  such  that   y(r)   is 

essentially  zero  for  targets  at  ranges  equal  to  or  greater  than 

one  or  two  track  spacings  away,  then  for  almost  all  values  of   i 

and   j,   p(i)   will  be  approximately  equal  to   p(j)   and  it  can 

rs/2 
be  assumed  that   p(det)  =  (N/Ns)      p(x;s)dx  where 

-s/2 

oo 

p(x;s)  =  l-exp[-  I     F(-js+x)]. 

j=-co 
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To  illustrate  an  application  of  the  above  result,  suppose 

Y(r)  =  kh/r3.   Then   F(x.)  =  2kh/wx2  .   If  it  is  assumed  that  the 

number  of  tracks  and  the  track  spacing  is  such  that   p(i)  =  p(j) 

for  almost  all  values  of   i   and   j ,   then 

m  °° 

F(x.)    =      I      F(-js+x)    =    (2kh/w)  (tt2/s2)csc2  (irx/s) 
j=-k  3  j=0 


so 


and 


p(x;s)    -    1  -  exp{- (2khir2/ws2)  esc2  (ttx/s)  }  . 


P(det)    =    2${(tt/2)  2    (W/s)}    -    1 

h 
where   W  =  2(2irkh/v)  '  is  the  searcher's  sweep  width  against  the 

target.   The  ratio  W/s   is  called  the  coverage  factor.   A  plot  of 

P(det)   is  given  as  a  function  of  W/s   in  the  National  Search 

and  Rescue  Manual. 

The  track  length  which  is  expended  in  an  unsuccessful  search 

is   nb    by  the  assumptions  of  the  parallel  sweep  model.    Since 

the  area  searched  is  nsb,   the  probability  of  detecting  a  target 

with  a  random  search  in  the  area  with   track  length  nb   is  given 

by   P(det)  =  1  -  exp{-W/s}.   This  value  could  be  compared  to  that 

given  by  the  parallel  sweep  model  in  order  to  obtain  a  relative 

measure  of  effectiveness  for  a  parallel  sweep  search. 
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X  .   The  Optimal  Allocation  of  Search  Effort  Problem 

The  problem  of  determining  an  optimal  allocation  of  effort 
for  a  search  is  discussed  next.  Here,  an  allocation  of  effort  is 
an  optimal  allocation  if  it  maximizes  the  probability  of  detecting 
the  target. 

The  class  of  searches  which  will  be  considered  are  those 
which  can  be  described  by  a  random  search  model  that  satisfies  the 
following  conditions:   The  target  is  fixed  in  a  region  of  area  A. 
The  region  consists  of  m   subregions  which  are  determined  by  the 
conditions  that  in  the   i —  subregion  the  sweep  width  is  a  con- 
stant w.   and  the  probability  density  function  is  a  constant 

p.  =  p./A.   where  A.   is  the  area  of  the  subregion  and   p.   is 
1    ri'  11  ^         *i 

the  target's  prior  probability  of  being  there.  (The  subregions 

are  assumed  to  be  numbered  so  that  W,  P,  >  W_P~  >  ...  >W  P  .)  If  a 

11    2  2         n  n 

target  is  detected  in  a  subregion,  the  searcher  must  be  in  the 

subregion  when  it  is  detected.   The  probability  of  detecting  the 

target  is  then 

n 
P(det)  =  I      {l-exp(-(p- )  >PiA. 
i=l 

where   (p.  =  W.A./A.   is  called  the  search  effort  density. 
i     ill 

Now  consider  the  problem  of  determining  the  solution  set 

cp.      which  maximizes   P(det)   subject  to  the  constraints 
^i  op 
^  n 

I   =   J  (A./W.)cp.  and  <p  .  £  0.   Note,   £   is  the  available  track 

i=l 
length.   This  is  a  nonlinear  programming  problem  for  which  the 

solution  set  is: 
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InW.p.  -(1/X)  7  (A./W.)ln  p.  +  £/A      i   in   ft 
11       ft  3      J  : 


(D  . 

yiop 


not  in   ft 


where   X  =  5!  A./W.   and   ft  =  {1,2,. ..,k}   where   k   is  determined 

ft   :   D 
by  the  condition  that  if   k  +  1   were  included  in   ft   then  <p,  ,  <  0 

Suppose  the  model  described  above  were  used  to  determine 

an  optimal  allocation  of  track  length  for  a  search.   The  result 

would  be  a  set  of  values   £«     =  (A./W.)<p-      which  maximize 

l  op      i   l  fiop 

the  probability  of  detecting  the  target.   Given  detection  will 
occur,  in  what  order  should  the  subregions  be  searched  in  order 
to  minimize  the  expected  time  until  the  target  is  detected?   The 
following  procedure  will  effectively  minimize  the  expected  track 
length  to  detection  and  gives  a  suitable  order  assuming  the 
searcher's  speed  remains  constant  so  that  minimizing  the  expected 
track  length  to  detection  is  equivalent  to  minimizing  the  expected 
time  to  detection. 

Divide   £   into  minimum  increments   A£   consistent  with 
the  random  search  model  and  then  allocate  the  increments  in  the 
following  order:   Allocate   Ail   to  the  Region  1.   This  will  be 
assumed  to  be  effectively  consistent  with  the  prescription  deter- 
mined by  the  formula  for   u> .      given   A£.   If  this  search  is 

J  Yiop   ^ 

unsuccessful,  determine  the  optimal  allocation  prescribed  for 
2A£.   This  will  effectively  indicate  whether  the  next  increment 
A£   should  be  assigned  to  Region  1  or  Region  2.   Continue  in  this 
fashion  until  detection  or  until   £   has  been  expended.   Note, 
this  is  an  approximate  procedure  in  the  sense  of  the  model. 
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For  a  discussion  of  the  above  procedure,  let   E(L.|l.££) 
symbolize  the  expected  track  length  to  detection  given  detection 
will  occur  by  the  time  track  length   £   has  been  used.   That  the 
above  search  procedure  will  effectively  minimize   E(l|l^£)   for 
an  optimal  allocation  of  track  length  is  suggested  by  the  follow- 
ing argument:   Divide   £   into  segments   A£   as  indicated  above, 
then  the  probability  that  the  target  is  detected  on  or  before  the 
i —  step  of  the  search   P(L^iA£)   is,  in  effect,  a  maximum, 
since  no  other  allocation  of  search  effort   iA£   will  give  a 
significantly  greater  value  for  any  value  of   i.   Hence, 
F  (iA£|L^£)  =  P (L^iA£ | L^£)   is  also  effectively  a  maximum  for  any 

Li 

value  of   i   since   P (L£iA£ |l*£)  =  P (L^iA£) /p (L^£)   and   P(L££) 

must  be  equal  to  its  maximum  value  for  any  alternate  procedure. 

n 
Therefore,   E(l|l^£)  =  I       [1-FT  (iA£ | L^£) ]   is  effectively  a 

i  =  l      L 
minimum  for  the  procedure. 

The  above  procedure  is  the  one  which  would  have  been  followed 
if  after  each  unsuccessful  search  one  calculated  new  prior  proba- 
bilities using  Bayes  procedure  and  then  did  the  search  prescribed 
for  a  track  length  A£  with  the  new  priors.  For  a  discussion  of 
this  point,  see  OEG  Report  No.  56.  An  application,  in  effect,  is 
discussed  in  Reference  4. 
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XI.   Target  Position  Probability  Distributions  Which  Change  in  Time 

Consider  the  following  model  in  which  a  target  moves  with 
constant  course  and  speed  in  a  plane  with  its  position  specified 
in  a  rectangular  coordinate  system:   The  target's  course  and  speed 
are  independent  random  variables   $   and  U  respectively 
whose  density  functions  £§($)      and   fn(u)   are  given.   The  joint 
density  function   fv  v(x,y;t)   of  the  target's  coordinates  is 

A  ,  I 

needed,  but  only   f    (x,y;0)   the  joint  density  function  for  time 

A  ,  I 

zero  is  given. 

For  such  cases,   fv  v(x,y;t)   can  be  found  in  principle  as 

A  ,  X 

follows:   First  it  is  convenient  to  use  polar  coordinates  and  with 
x  =  p  cos  a   and   y  =  p  sin  a   and  to  set 

f(p,a;t)  =  fx  y(p  cos  a,p  sin  a;t)  . 

Note,  to  first  order,   P{target  is  in   AA  at   t}  =  f (p,a;t)AA. 


x 


Figure  3.   Problem  geometry 
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Referring  to  Figure  3, 


f  (  p  ,  a ;  t )  A  A  = 


2tt 


f  [r(u,<|>)  ,0(u,<|>)  ;0]AA,f$(<|))fu(u)d«J)du 
0   0 
where  the  functions   r(u,<j>)   and   6(u,<j>)   are  such  that  if  the 

target's  polar  coordinates  were   r(u,<}))   and   0(u,<}>)   at  time   0, 

then  with  the  speed  and  course   u   and   <J>   they  would  be 

(p,a)   at   t.   The  functions  can  be  expressed  implicitly  by 

r2  =  p2  +  (ut)  2  -  2p  ut  cos  (a-<j>)   and   r  cos  9  =  p  cos  a  -  utcoscj). 


The  integration  is,  in  effect,  a  sum  over  all  possible  pairs  of 
values  of   u   and  $ .   Each  pair  determines  a  vector   ut  which 
translates  target  positions  at  time   0   to  target  positions  at 
time   t   such  that   AA   is   AA'   translated  without  rotation  or 
distortion.   So,   AA1  =  AA,   independent  of   u   and   <J> .   Therefore, 


f(p,a;t)  = 


2tt 


f  [r(u,<f>)  ,9(u,<|>)  ;0]f$((J))fu(u)d4)du 


0   0 
As  an  example,  suppose 


(x2+y2) 


fxyx,y;0)  =  2^-  e 

with   f*(4>)  =  -i-  ,   0  <  d>  £  2tt,   but  with   U  =  u   known.   In  this 
case , 

_  r2(u,(f>) 
f  [r(u,<t>)  ,6(u,<j));0]  =  — ^  e 


and,  using   r2  =  p   +  (ut)  z  -  2p  ut  cos  (a-<$>)  , 
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f(p,a;t)    = 


•2tt    -   ^p-    [p2+(ut)  2-2p  ut  cos(a-<J>)  ] 


2ttcj 


d* 
2tt 


2ira 


t  e 


2a 


*-    [p2+(ut)2] 


_1_ 

2tt 


n„    put  ,        ,  . 

f2ir   i^-r-  cos  (a-(J>) 


d<J> 


-r    [p2+(ut)2] 


2tra 


■r  e 


2a'    tr-    ■  v~w/    J         riPut^ 


where   Jn   is  the  zero  order  Bessel  function  of  the  first  kind 
and    i  =  /-l.      To  obtain   f    (x,y;t),   substitute   P2  =  x2  +  y2 

A,  I 

in  the  above  result.   In  this  case,  since   fx  y(x,y;t)   is  symmet- 
ric about  the  origin,  the  marginal  distribution  of  the  target's 
range   R   from  the  origin  is  of  interest.   In  general, 


f   Q(P,a;t)  =  pf  y(p  cos  a,p  sin  a;t)  =  pf(p,a;t) 


and 


fR(P;t)  = 


27T 


fR  0  (p,a;t)da 


so,  in  this  case, 


fR(P;t)  =  £-e 


2^  tP2  +  (ut)2] 


The  distribution  function   F  (p;t)   for  this  case  is  given  in 
Reference  5. 
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